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Abstract 

The basics of the Wigner formulation of Quantum-Mechanics and few related inter- 
pretational issues are presented in a simple language. This formulation has extensive 
applications in Quantum Optics and in Mixed Quantum-Classical formulations. 

1 Introduction 

In classical mechanics, a system is completely specified by the position and momentum of its 
particles. It is common to build a 6N dimensional space, where N is the number of particle 
constituting the system. This space is called the phase space and contains an axe for each of 
the coordinate and momentum of the system. The system is completely specified by a point 
in this space. Usually, when we have a lack of knowledge regarding the state of the system, 
we try to find a probability distribution for the system in its phase space. This is the case in 
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classical statistical mechanics, where f(p, q) (p is a 3 AT dimensional vector representing the 
momentum of all particles, and q is a 3 AT dimensional vector representing the coordinate of 
all particles) is the probability of finding the system with coordinate between q and q + dq 
and momenta between p and p + dp. The predictions of quantum mechanics are similar to 
that of classical statistical mechanics in the sense that they are statistical in nature. Thus, 
the fundamental question arise whether quantum processes can be described as an average 
over uniquely determined processes or not. And the observability of these processes if the 
answer is yes [T]. Naturally the place to look for such processes is the phase space and if 
the system is undergoing a process in the phase space it must be represented by a point 
in the phase space. Regardless of the tough problem mentioned, it seems natural for a 
scientist to seek a function similar to distribution functions in quantum mechanics. But, in 
the framework of Orthodox and Copenhagen interpretations of quantum mechanics this is 
impossible. 



2 Wigner distribution function 

The first example of such a function in quantum mechanics was suggested by Wigner [2] . He 
mentioned the canonical ensemble probability distribution in classical statistical mechanics 
for a system having momenta between p and p + dp and coordinate between q and q + dq, 
which is e~P e , where f3 is the reciprocal of the temperate, T, and e is the sum of the 
kinetic and the potential energy. He mentioned that in quantum mechanics, we cannot 
simultaneously have momentum and position so we cannot have such an expression. But 
even if we consider the coordinate alone, where the classical expression for probability is 
e~@ v , where V is the potential energy of the system, the classical expression is not valid for 
quantum systems, because when (3 — > oo there is no reason for that expression to be equal to 
\ipo(xi, x n )\ 2 (the ground state wave function is not even always known). Although, the 
statistical mechanics of quantum systems is given by the von Newman formula, i.e., < Q >= 
Tr(Qe~" ), where Q is the operator corresponding to the quantity under consideration, H 
is the Hamiltonian of the system and <> denote the expectation value. Since, it was not 
easy to use the von Newman formula for evaluating the expectation values, Wigner suggests 
to build the following expression 

P(q,p) = (^-) n / dytftq - y/2)^(q + y/2)*e^\ (1) 

and call it the probability function, here n is the dimension of the space. Thorough out 
this paper consider the limits of the integrals from — oo to oo, unless otherwise is explicitly 
stated. Unfortunately, Wigner never mentioned how he come up with this recipe, he just 
mention: "This expression was found by L. Szilard and the present author some years 
ago for another purpose." By introducing the inverse Fourier transform of tp, i.e., tp(q) = 
2 / dpe^ pq ip{p) m the above relation we could get 

J dy J dp' J dp"^{p'y^(p")e^ [ - p ' il]+v/2)+p " {l] - y/2)+py \ (2) 

by performing the integral over y we could get (2wh) n S(p — p ^ p ), then we can perform 
the integral over p" and perform the change of variables p — p' — > —y/2to get 

P(q,p) = (t^-Y 1 I dy^(p - y/2M P + y/2)*e^ h . (3) 
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This relation is completely equivalent with the relation |T]) and shows the symmetry of the 
Wigner functions with respect to q and p. The phase space function corresponding to an 
operator A is defined thorough 

A(q,p) = f dye^ h < q - V -\A\q + V - > , (4) 



3 Proposals for getting the Wigner function 

Stenholm presents a derivation for the Wigner function [3] . All the information extractable 
from the quantum theory is contained in the matrix elements 

< xi\p\x 2 >= ip(xi)i>(x 2 )* ■ (5) 
We can bring the density matrix into momentum representation and write 

< Pi\p\P2 >= I f dx 1 dx 2 exp[-i(p 1 x 1 -p 2 x 2 )/h] < Xi\p\x 2 > ■ (6) 



Similar to a two body problem in mechanics, we can define new variables as R — Xl ~^ X2 
and r = x\ — x 2 , and a similar change of variables in the momentum representation, i.e., 
p _ pi+P2 an( j p — p l ~ p 2 . It is simple to show that 

P\X\ - p 2 x 2 = Pr +pR. (7) 
By substituting ([7]) into ([6]) and changing the variables, we could get 

<P + \\p\P~\ >=^[ J drdRexp[-i(Pr+pR)/h]<R+ r -\p\R- r -> . (8) 

The above relation is just the Fourier transform of p{R 1 r) 1 where by analogy to the two 
particle problem, we can call R the center of mass coordinate, and r the relative coordinate. 
Because, we are interested to get a function containing both momentum and coordinate, 
we could either drop the Fourier transformation on relative coordinate to get the Wigner 
function, or drop the Fourier transformation on the center of mass coordinate to get the 
Shirley 4 function. 

Groot has presented another equivalent method for deriving the Wigner function 
By inserting unity operators, we get 

A = j dp'dp"dq'dq"\q" >< q"\p" >< p"\A\p' >< p'\q' >< q'\. (9) 



Then, we can introduce the new variables p' = p — u/2, p" = p + u/2, q' = q — v/2 and 
q" = q + v/2, where the Jacobian is equal to one, and use the relation < q\p >= h^ n / 2 ei p ' q 
to get 

A = h~ n J dpdqdudv\q + v/2 > e i(p+«/2)(<?+«/2) <p + u /2\A\p - u/2 > 

(10) 

e -fr(p-u/2){q-v/2) <g _ v / 2 |. 

This relation simplifies to 

A=hT n J dpdqdudv\q + v/2 >< p + u/2\A\p - u/2 X q - v/2\ex( qu+pv l (11) 
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By defining the A dependent function 

a(p, q)= J du<p + u/2\A\p - u/2 > e% qu , (12) 
and the A independent operator 

A(p,g) = Jdv\q + v/2 >< q-v/2\e% pv . (13) 

We have 

A = h~ n J dpdqa(p,q)A(p,q). (14) 

It is clear that a(p, q) is the Wigner function corresponding to the operator A. This is a 
natural way one can lead to the definition of the Wigner function. 

4 Weyl operator 

Before this work by Wigner, Weyl [6] had proposed a method to construct an operator A 
corresponding to the phase space function A(q,p). First we define 

a(a,r) = C^A J dq J dpe~^ + ^' h A(q,p) (15) 

and then, 

A(q,p)= j da j dTa(cr,T)e l ^ +T P^ h (16) 



Wigner's recipe is exactly the inverse of the Weyl's. If this is a suitable correspondence 
between A(p, q) and A, so we must be able to get the correct expectation value for A by 
use of A(p, q), i.e., 

< ^\A\4> >= [dq [ dpP(q,p)A(q,p). (17) 



Before proving this equality, I should mention a lemma. 

lemma 1: By using the Baker- Hausdorff lemma, we can prove that 

e A+ * = e A e*e^ A > B \ (18) 

which yields to 

e i( CT, 3+ T p) = e jt<rq e i r P e i<rT / 2h , (19) 
By substituting A(p, q) from (fT5|) and A from (fj"6|) into (fT7|) . we get 

J da J dra(a,T) < i}j\e l ^ +T ^/ n \^ > 

= JdaJdrJdqJ dpP{q,p)e^ aq+T P^ h a{a 1 T) 1 
which easily simplifies to 

< ^\e l ^ +T ^' h \^ >= Jdqf dpP(q,p)e i ^i+ T P^ h 
= {2irh)- n JdyJdqJ dp^(q + y)*i/j(q - y) e i ( 2 PV+^+r P )/h^ 



(20) 



(21) 
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The integral over p gives {2irh) n 5{2y + r), which allow us to perform the integral over y in 
order to get for the right hand side 



dqip(q + r/2)*V(g - t /2)e t(aq)/h . (22) 
According to the lemma 1, the left hand side is 

Because p is the generator of translation (j2"3"|) is equal to 

dxe^^ w+aT ^^( X ytp{x + T). (24) 
By imposing the change of variable x — > q — r/2, we get the relation (f22|) . so Q.E.D. 

5 Properties of the Wigner distribution 

A number of properties have been mentioned for this function |7] 

(i) Since P{q 1 p) should be real, it should be corresponding to a Hcrmitian operator, i.e., 

P(q,p)=<^\M{q,p)\^>, (25) 



where M — Aft, i.e., Hermitian. 
(ii) 

a ipy/h 

(26) 

= / dyd{y) <q- y/2\p\q + y/2 >= \iP{q)\ 2 =< q\p\q > . 



J d P P(q,p) = (5^)» Jdpjdy<q- y/2\p\q + y/2 > e l - 



dqP(q,p) = \ij(p)\ 2 =<p\p\p> (27) 

JdqJ dpP(q,p) = Tr{p) = 1 (28) 

Derivation of the second and the third one are similar to that of the first one. 

(iii) Translation of P(q,p) in the momentum and coordinate spaces occur in accordance 
with the translation of the wave function, i.e., if ij)(q) — * ip(q + a ) then P(q,p). — * P(q+a,p), 
and if ijj(q) — > e lp q / h ip(q) then P(q,p) — > P(q,p — p') (iv) P{q,p) should change the same 
way as ip m space reflections and time inversions , i.e., if ip(q) — > ip(—q), then P(q,p) — > 
P(-Q, ~P) and, if ijj(q) -> ^(<?)* then P(g,p) P(<?, -p) 

(v) When the third and all higher order derivatives of the potential are zero we get the 
classical equations of motion (the Liouville equation). (This will be shown in the section 
Dynamics.) 

(vi) 

| < i>(q)\4>(q) > | 2 = 2nh JdqJ dpP i! (q,p)P (t> (q,p) (29) 

(vii) 

J dqj dpA(q,p)B{q,p) = 27r?iTr(AB), (30) 

where A(q,p) is the classical function corresponding to the quantum operator A. Using 
the property (ii) it can be easily shown that if h(q,p) = f(q) + g(p) then we can get the 
expectation value of h by J J dpdqP(q,p)[f + g] 
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6 The Product of two Operators 



Groenewold in a fundamental work presented some foundational issues of quantum mechan- 
ics. He depicts the physical properties corresponding to the quantum mechanical operators 
A and B with a and b. He used the von Newman's assumptions, i.e., (I) if a corresponds 
to A and b corresponds to B then a + b corresponds to A + B, and (II) if a corresponds to 
A then f(a) corresponds to f(A). He shows that such symbols constitute two isomorphic 
groups. Thus, if A and B do not commute then a and b should not commute. It can be 
shown that by assuming a and b as commuting observables we get into contradiction with 
assumptions (I) and (II). Therefore, a quantum system can not possess two physical prop- 
erties corresponding to two non-commuting operators, and there is no reason to introduce 
different notation for operator and physical property. In that paper, he shows 

AB = F - F(q,p) = A(q,p)e^ A /^B(q,p) = B(q,p)e^ hA /^ A(q,p), (31) 

where < ^ < 

d d 3d 

dp dq dq dp ' 

is the negative of the Poisson bracket. Note that there is a dot product between the 
differentiation toward right and differentiation toward left. By taking the matrix elements 
of (fPBj) . we get 

< q"\A\q' >= J da J dra(a,T) < q"\e l{all+T P )/h \q' > . (33) 

and by using the lemma 1, we can get 

< q"\A\q' >= J dadTa(a,T)e' aT ^ n e ia ^'- T ^ h 5{q' - 
= jdaa{a, q' - g ") e M<?'+?")/2&. 



q") 



(34) 



Now, we have 

F(q,p) = J dz^ z / h <q-l\AB\q+Z > 
= J dzdq'e lpz / h < q- §\A\q' >< q'\B\q+^ > 
J dzdq'dada'e^/^^'+i-^e^/^'^'+i+^aia, q' - q + %)a'(a', q-q' + ^)e zpz/h . 

(35) 

By defining the new variables r = q' — q + | and r' = q — q' + # , we would get 
F(q,p) = f dTdT'dada'a(a,T)e^ h '> (aq+T P^e^ 2n ^ a ' T - aT ^ 

(36) 

Xe (i/h)^' q +r' P ) a f^yy 

Now, consider the exponential between the other two exponentials and Taylor expand it. 
Consider the second term while forget about all constants, i.e., <t't — or' , it is easy to see 
that by differentiation of the exponential on the right with respect to p and the exponential 
on the left with respect to q, we can get er'r. We can get err' by differentiation of the 
exponential on the right with respect to q and the exponential on the left with respect to p. 
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Therefore, replacement of (i/h)(a'T — <jt') by (HA/2i) makes no difference up to the second 
term in the Taylor expansion, by more elaboration you can show that this is also true for 
the higher order terms. After the mentioned replacement A(p, q) and B(p, q) (inverse of 
(p"5"j) ) will appear in (|3U)) , and we will get the first equality in (|3"Tj) . If we change the place 
of the first two and the last two terms in (|36[) . we can repeat the preceding discussion by 
interchanging the differentiation with respect to p by the differentiation with respect to q 
and vice versa. Thus, we can easily get the second equality in (f3T|) . 

Another way of writing the product of two operators are Bopp operators which are 
defined as [5] 

^ h d „ h d .„„, 

Q = 1~^^ P = P+^7r- ( 37 ) 



2i dp ' 2i dq 

dq U dp J 



By taking a test function / and a little elaboration you can show that [aq + rp, ~~ <J ~§^\ = 
0. This equality yields 



"*> {1 w« = «* ( ~'"« 1< '*-'* , • <»> 

If we multiply both sides by e^^ CT q+T p \ Taylor expand the middle term on the right hand 
side, and operate it on the exponential on its right, then every will be replaced by r' and 
every ^ will be replaced by & ' . Then, we will have the Taylor expansion of an exponential 
function in the middle, which can be gathered and give the final relation 

(39) 

= e w( a i+ T P) e 2k( rcr ' ' - aT ') e i( <T ' 'i+ T ' p) . 

On the right hand side of (f39|) , we have all the exponential terms we had on the right hand 
side of (|36p . By replacing them and using the notation introduced in (|37[) . we get 

F{q,p) = J dTdr'dada'a(a, T)ei^ +TP ^ei^'i +T 'P^a'{(j' , t'). (40) 

Now, we can define 

A{Q,P) = J dTdaa(<j 7 T)e*(' 7Q+TP l (41) 
Therefore, we can express F(p, q) as 

F(p,q)=A(Q,P)B(p,q). (42) 

Similarly it can be shown that 

F{p 1 q) = B{Q\P*)A{p,q) 1 (43) 

where 

n* , h 9 v* % d 1AA\ 

Q =q+ 2ldp-> P =p -2id q - (44) 

We know that the Wigner function is the function associated with {^pg) n p, and that the 
equation of motion for p is 

ihdp/dt=[H,p}. (45) 
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Using the product rule, we can transform (T4"5|) to 

ikdP/dt = H(q,p)e hA / 2l P{q,p) - P(q,p)e hA ^H(q,p). 



(46) 



The first term of the Taylor expansion is HP — PH, which is equal to zero. The second 
term in the Taylor expansion of the first and the second term of (|4"6"|) are just negative of 
each other so they build up to — ^f"§f ]• For the third term we need 



^2 j9 9 d a ^.d d d d ^ 
dp dq dq dp dp dq dq dp 



(47) 



By inserting two test functions / and g, respectively, in the left and right hand side of the 
expression in (|47p . we can show that 



A 



P P 



P P 



P P 



dp 2 dq 2 dqdp dpdq dq 2 dp 2 
In order to generalize the expression for A 2 to higher dimensions, we can write 



A 



d d 



d d 



dpi dqi dqi dp t 



d d 



d d 



dpj dqj dqj dp. 



in order to get 



A2 = E 



d 2 



P 



d 2 



P 



d 2 



P 



dpidpj dqtdqj dq t dpj dpidqj dqtdqj dpidpj 



(48) 



(49) 



(50) 



Because of the symmetry of (148I) . the third term in (|45|) is zero. Again, by inserting the 
test functions / and g, we can evaluate 



A 3 =A 2 A = f^ ? _2-SL_2_ + & $)(JLJL 

V Q v '2 Q a 'i dqdp dpdq dq' 2 dp' 2 ' ^ dp dq dq 



d d 



P P _ 3 P P 

dp 3 dq 3 dqdp' 2 dpdq' 2 



d 3 d 



7$ 



' Qq'2Qp dp' 2 dq dq 3 



P P 



and 



A 4 = A 3 A 



P P 

dp 3 dq 3 



3 b 3 P +3 d 3 ,7 

dqdp' 2 dpdq' 2 dq' 2 



P P 

dp 2 dq dq 3 dp 3 



' 4 P _ 4 P P | q d 4 

dp^ dq A dqdp 3 dpdq 3 dq 2 dp 2 dp 2 dq 2 



d 4 d 4 



P 



Z73 



4 d 4 & 



/ d 8 _ d d \ 

^ dp dq dq dp > 



d 4 d 4 



dq A dp dp^dq dq A dp 4 

By continuing in this manner we can show that 

hdP/dt = -2H(q,p)sin(hA/2)P(q,p). 



(51) 



(52) 



(53) 



7 Proof of the impossibility of a positive phase space 
probability distribution 

Wigner was aware that this probability function gets negative values (unless the world 
was just made up of Gaussian wave packets). Thus, he emphasized that this is just a 
calculational tool, not a real probability distribution in phase space. Latter, he uses the 
fact that for a mixed state P(q,p) = ^2u>iPi(q,p), where Wi is the probability of the i'th 
pure state, and Pi(q,p) is the Wigner function for the i'th pure state. Then, he shows that 
by imposing the conditions (i) and (ii) (see section l[5])) it is impossible to build an always 
positive distribution function [5]. He used ip(q) = aipi(q) + bip 2 (q), where is zero outside 
I\ and ip2 is zero outside of I 2 . Consider I\ and 1% to be two non-overlapping intervals over 
the space of coordinate. Now we have 

P ab (q,p) = \a\ 2 Pi + a*bP 12 + ab* P 21 + \b\ 2 P 2 . (54) 

If q is outside of I 1: Pi is zero for such a q and the only way to have a positive value for 
P a b{<liP) for every a and b is to have -P12 (<?,£>) = P2\{<l,p) = 0. The same reasoning can be 
given for the q outside of I 2 . Therefore, every where, we have 

P ab {q,p) = M 2 Pi + \b\ 2 P 2 . (55) 

This means that P a b is independent of the complex phase of a/b. Consider the Fourier 
transform of ipi and ip2 to be <p\ (p) and <f> 2 (p) ■ By removing P a b from equation (|54[) and 
equation (|55[) . then integrating both sides of the resultant equation with respect to q and 
using the mentioned Fourier inverses, we get 

\a\ 2 JP 1 (q,p)dq+\b\ 2 JP 2 (q,p)dq 

(56) 

- |o| a |0i(p)| 2 + \b\ 2 \Mp)\ 2 + 2Re[ab*Mp)MPn 
For this relation to be valid for all a and b, we must have 

Mp)Mp)* = 0- (57) 

On the other hand, 4>\ and 4> 2 are Fourier transforms of confined functions; thus, they 
cannot vanish on a finite interval. This is a contradiction, and QED. Because, it seems 
possible to break down every normalizable wave function into such a linear combination, 
thus this proof excludes the possibility of having a phase space distribution for a quantum 
state. Wigner [9], also, showed that by imposing the conditions (i)-(v) |T]) is unique. While, 
O'connell and Wigner [10] show that by imposing conditions (i)-(iv) and (vi) is the only 
possible distribution. 



8 Dynamics of the Wigner function 



If we want to express quantum mechanics in terms of the Wigner function we must derive 
Wigner functions equation of motion. This will be done with the aid of the Schroedinger 
equation, i.e., 

. dt[>(t) 



ih- 



dt 



2m dq 2 



V(q,t) 



(58) 
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By conjugate transposing both sides of the Schroedinger equation, we get 

h 2 d 2 



in- 



ch 



2m oq 



Decomposing the time dependence of P into two parts, we have 

w - (wT I d y f 2 *^^ -v) + *(« + y)* 9 -^ 1 



a 2ipy/K 



(59) 



(60) 



9k P i d v P 
dt T dt ■ 



In the last expression of (|60|) . the first part arises from the kinetic part of the Hamiltonian 
and the second part arises from its potential part. By substituting (|58|) and ([59]) in (|60|) 
and considering the n to be equal to one, we can get 



dt V 27rm 7 ' 



5 — Hq-v) -i'w + v) — ^ — 



„2ipy/h 



(61) 



where we have replaced d 2 /dq 2 by d 2 /dy 2 . Integration by parts, because ip vanishes at -co 
and oo, yeilds 



— = ( ~ P ) / 



^ + y) %(g-»)-^ + yr^ (? - y) 



a 2ipy/h 



By going back to d/dq, we obtain 



dfc-P _ p dP(q,p) 
dt vci dq 



(62) 



(63) 



which is identical to the corresponding term inthe classical Liouville equation. Also, we 
have 



%f = / WV>)(? + vTHq -y)~ 1>(q + v) - y)]e**»/* 

= / dv[V{q + V)- V(q - y))^{q + y)*^(q - y)e 2 ^/ h . 

By Taylor expanding V, we get 



(64) 



A=0 



A! <v 



Therefore, we have 



a,P 2i 



dt TTh 2 I 4 ^ 



y x d V 
AT TV 



VX<z + y)*VK<? - y)e 2w/r \ 



(65) 



(66) 



where the sum is over the odd positive integers A, since the even terms resulting from 
V(q + y) and those resulting from V(q — y) cancel each other. Because by differentiating 
the exponential term with respect to p, we get a y multiplier, y x can be replaced with 
[(K/2i)(d/dp)] x to get 



dq x dp x 



(67) 
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For the sake of simplicity, equations (|65[) - (pIT|) are written for the one dimensional case 
(n = 1). In order to generalize them to higher dimensions A! should be replaced by IFA^!, 
any thing to power A with the same thing to power J2i ^i; 9q x with Uidqf, and dp x with 
Tlidpj' . But, remember Ai's take positive integers that yield an odd positive integer for 
J2i ^i- We can, also, write ^f- in the form 



d v P 
dt 



djP(q,P + j)J(q,j), (68) 



while 

J(QJ) = JWy^K I dy[V(q + y) - V(q - y)\e^l n 

(69) 

= T^rn I d y[ V (<l + V)~ V (l ~ y)]sin(2jy/h) 

has interpreted as the probability of a jump in momentum by an amount j, if the position 
is q. We can go from the first to the second equality in (p9")h because e %x = cos(x) + isin(x) 
and the function in the square brackets is an odd function so when it is multiplied by an 
even function, cos(2jy/h) and integrated over the whole space, it will give zero. 
Now, we are able to get the equation of motion as 

dp _ op v dE^v (ft/^E*-- 1 cg>p 

8t 2^ m t dq % + ^ Utdq* U\\ Hdp^ ' 1 ' 

Now consider the case where the potential has no third or higher order derivative; then, 
evidently ([M]) is the classical Liouville equation. For a system consisting of a bunch of 
harmonic oscillators and free particles, surprisingly we can solve the easier classical equations 
of motion and get the exact quantum result! 



9 Attempts for giving a probability distribution inter- 
pretation to the Wigner function 

Some people have argued this function as a valid probability distribution, and some others 
argue it as a valid probability distribution just for some situations. Stcnholm has argued 
that, we can "obtain verifiable predictions" only by using "suitable test bodies." He show 
that always positive probabilities come out of the Wigner distribution when these arguments 
are implemented. In relativistic quantum mechanics }llj . even if we are working in the 
position representation, the choice of the position observable is not at all trivial. Therefore, 
probably in non-relativistic cases it is just the absence of mathematical complications which 
make us to believe that we can make a classical interpretation of position and momentum. 

The way to measure the momentum and coordinate of a particle is to let it interact with 
another body which we usually let approach it's classical limit. This second body is a test 
particle, and we are actually performing a scattering experiment. In this scattering exper- 
iment, test particle transfers the desired information out of the interaction region. "Only 
probability distributions observable in this manner can be given a physical interpretation" 
[3]- In cases where the test particle is carrying both coordinate and momenta information, 
restrictions due to Heisenberg uncertainty principle should be taken into account. Stenholm 
emphasized that, in order to get the most precise results for both momentum and coordi- 
nate, we must use a test particle which is in a state of minimum uncertainty. At the end, 
the test body is bringing out some information which at best allow us to confine our system 
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to a region of phase space satisfying the relation ApAq < | and no more precise detail is 
achievable. The minimum uncertainty wave packet can be determined uniquely as 

/ (x-<x>) 2 i<P>x \ 
ip (x) = Cexp I — 2 + I , (71) 

where <> denotes the expectation value, and we have uncertainties Ax = b and Ap ~ h/2b. 
The Wigner function for this minimum uncertainty wave packet is 



, (R- < x >) 2 46 2 (P- < P >f , 
W (R, P) = Aexp ( ^-J- [ —^- 2 ) • ( 72 i 



Stenholm argued that "a Wigner function W(P, R) is not directly observable but has to be 
convoluted with the function describing the test particle, which smears it, at least, by the 
amount implied by the function" (j?2"|) . This convolution leads to 

P(tt, q) = A J J exp (~^# - 2b ' 2{P h ^ )2 ) W(R, P)dRdP 

(73) 

= {£r)ffj drdPdRe^ Pr ^e^ R ~^ 2 / 2b2 e- 2b ^ p -^ 2 / n2 i;(R + §)V(i? - §)*. 

Here (q, ir) are a couple of position and momentum coordinate, and Stenholm hope them to 
give meaningful phase space interpretation. A is just a normalization constant. By carrying 
out the integral over P, we can show that P(q, 7r) > 0. The immediate criticism to this 
approach is that not all test bodies are minimum uncertainty wave packets. Thus, Stenholm 
emphasize that in order for the test body to exhibit a nearly classical behavior it must be 
smooth and every wave function smoother than the minimum uncertainty wave packet will 
fulfill the positiveness requirement. All in all, we have to calculate W(P,R;t) up to the 
moment of measurement, and then smooth it to obtain P(tt, q; t). Since, P(ir, q; t) depends 
on the test particle prepared by the observer it has no dynamics. The idea of smoothing the 
Wigner function with a Gaussian function was first introduced by Husimi p~2] . He get the 
positive distribution which is now called the Husimi distribution. Husimi didn't interpret 
it as a phase space distribution because it doesn't poses the property (ii). There is also 
a bunch of work for interpreting Quantum Optics based on Wigner function. Including 
Marshall and Santos [13] . who argue that, there is just a subset of states in the Hilbert 
space which can be generated in the laboratory, and those states "May be represented by 
a positive Wigner distribution." They have also claimed that [2], the experiments which 
are exhibiting non-classical behavior of light can be interpreted just by assuming light 
as an electromagnetic wave in accordance with the Maxwell's equations of motion. Also, 
Holland et al [T5] published on the " Relativistic generalization of the Wigner function and 
its interpretation in the causal stochastic formulation of quantum mechanics." There are 
some good reviews on the mathematical properties and applications of the Wigner function, 
e.g., m, M and [17]. 
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We can argue that, any experiment which is designed to measure both momentum and co- 
ordinate will get average information about a region of phase space which is large enough to 
give a positive value. The Wigner function give a positive average over this region and thus 
interpreting the Wigner function as a probability distribution function is experimentally ad- 
equate. But does it make sense to consider a probability function which has no meaning on 
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a point but is representing a physical reality when averaged over any large enough interval! 
Thats weird but is it more weird than a particle owning either momentum or position and 
not both of them at the same time? The classical formalism of quantum mechanics doesn't 
provide any prediction regarding a simultaneous measurement of position and momentum. 
If it was possible to perform such an experiment, for example, as suggested, thorough send- 
ing out a test particle to interact with system (quantum particle) and be able to measure its 
position and momentum before and after the interaction, because it is approaching classical 
behavior, then we could perform an experimental test for this interpretation. But it seems 
that no semi-classical test particle can be used for a measurement on a quantum system 
without enormously changing its state. 
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